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Fermi gases with generalized Rashba spin orbit coupling induced by a synthetic gauge field have 
the potential of realizing many interesting states such as rashbon condensates and topological phases. 
Here we develop a fluctuation theory of such systems and demonstrate that beyond-Gaussian effects 
are essential to capture the physics of such systems. We obtain their phase diagram by constructing 
an approximate non-Gaussian theory. We conclusively establish that spin-orbit coupling can enhance 
the exponentially small transition temperature (T^) of a weakly attracting superfluid to the order 
of Fermi temperature, paving a pathway towards high Tc superfluids. 
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Construction and study of model Hamiltonians with 
quantum gases has opened up the possibility of not only 
addressing long standing questions [THS] but also creat¬ 
ing systems that are not conventional. The recent ad¬ 
vances in synthetic gauge fields [7ffT^ have provided new 
impetus, motivating studies of interacting bosons and 
fermions in their presence that are of interest to a wide 
array of physicists (see review Ref. IT5)) . 

A uniform non-Abelian gauge field results in a gener¬ 
alized Rashba spin orbit coupling (RSOC). Interacting 
fermions with RSOC have many interesting and novel 
features, [iMni while additional Zeeman fields can help 
realize topological states. [TBl - BD] Even in the presence of 
weak attractive interactions, a crossover from a BCS type 
superfluid to a rashbon-BEC can be achieved by increas¬ 
ing the strength of RSOC. mus] Rashbon-BEC is a con¬ 
densate of rashbons - a bosonic bound state of a fermion 
pair in presence of large RSOC - whose mass (between 
2-2.5 fermion mass) determines its transition tempera¬ 
ture which is of the order of the Fermi temperature. In¬ 
deed early studies[Tni[lI] suggest this enhancement of the 
transition temperatures of weakly attractive systems by 
means of RSOC. In addition to this, the rashbon conden¬ 
sate/gas has several uncommon and unique traits - for 
example, unlike the case of the boson-boson interaction 
in the usual BEC,[22] the rashbon-rashbon interaction is 
independent of the interaetions between the constituent 
fermions.^^ The hunt for the ‘best possible Tf in these 
systems, along with the novel physics just noted, mo¬ 
tivates this study to take the vital, if challenging, step 
- construction of a finite temperature theory including 
fluctuations beyond the mean-field. 

In this report, we develop, for the first time, a fluctua¬ 
tion theory of the normal state of an interacting Fermi gas 
with RSOC. We show that the Gaussian theory, which 
provides an excellent qualitative description of the BCS- 
BEC crossover in systems without RSOC P^IBS] is woe¬ 
fully inadequate to describe the normal state of systems 
with RSOC even at a qualitative level. We develop an ap¬ 
proximate theory, including the crucial beyond-Gaussian 
effects, and use it to obtain the phase diagram of inter¬ 


acting Rashba Fermi gases. Novel results include a clear 
demonstration of the enhancement of the superfluid tran¬ 
sition temperature (Tg) in weakly attracting system from 
an exponentially small value to that of the order of Fermi 
temperature (Tp). We also show that in the regime of 
weak interactions, the superfluid transition temperature 
is a non-monotonic function of RSOC. 

Formulation: Choosing units where the fermion mass 
m and the Planck’s constant h are unity, the kinetic en¬ 
ergy of a spin-orbit coupled fermion with momentum k 

in three spatial dimensions is efco- = ^— a|fcA| H— f-- 
Here a{= ±1) is the helicity, k\ = Xxkx^x + XykyBy + 
Xzkz^zi where A = {Xx,Xy,Xz) = X\ describes a ‘vec¬ 
tor’ in the gauge field configuration space, and Xm = 
Max(Aa;, Ay, Ax). Non-Abelian gauge fields(equivalently 
RSOC) with high symmetry, such as the spherical gauge 
field with Xx = Xy = Xz = A/a/ 3, are of particular in¬ 
terest. We refer to the absence of gauge fields/RSOC 
(A = 0) as “free vacuum”. 

A finite density po of fermions determines a charac¬ 
teristic momentum scale kp defined by po = kp/iir'^ 
and an associated energy/temperature scale Ep = Tp = 
kp/2. The singlet interaction (bare strength v) between 
the fermions is characterized by a scattering length Og. 
Physics at temperature T and chemical potential p with 
volume V is studied using functional integral methods. 
After introducing pairing fields r]{q),q = (iqi,q) (iqi- 
Bose-Matsubara frequency, q ~ wave vector) and inte¬ 
grating out the fermions, the action 

S[q, F]=- lndet[-G-i] - ^ E, E(g) oid)- (1) 

is obtained, where F is a source field and G is the Greens 
function (functional of rj and F).[2^ To study the normal 
state physics, we expand the exact action Eq. Q about 
the saddle point where q{q) = 0, 

5 « _ lndet[-Go-i] - i E, V*{q) vid) + E, lHd)L{d)lid) 

+ Eq1.g2.93.94 'y*idih*id 2 )K{qi,q 2 -, qs, 94 ) 7 ( 93 ) 7 ( 94 ) ( 2 ) 

up to quartic order (7 = rj+F)-, Go is the non-interacting 
Greens function. The quantities L and K are derivatives 
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of the action (Eq. Q) to appropriate order in rj. They 
are constrained by conservation laws, e.g., the arguments 
of K have to satisfy momentum conservation. 

Gaussian Fluctuation Theory: Retaining only the 
first three terms in Eq. ([^ produces the Gaussian fluc¬ 
tuation theory [241 - 126] . quadratic in 77 . Upon integration 
of the rj helds, we obtain 


Sg[F] = -lndet[-Go^]-bX; 5 lnM(q) F*{q)x{q)F{q), 


(3) 

where M(g) = L{q)-^ = Hq) + vJ2k Hq) = 






iqi-i 


(2+fc)o 


with np{x) = -I- 1) denoting the Fermi function, 

Aap is the amplitude of the singlet in the two particle- 
state with momenta q/2 ± k and helicities a and /3, 
and ^ = e — fi. The analysis also produces the pairing 


susceptibility x(g) = L{q) ~ 1^? whose divergence 

from the positive side up on the reduction of tempera¬ 
ture indicates a pairing instability. The first such diver¬ 
gence of x( 0 , q) occurs at q = 0 (as we have verified), 
i.e., the system is most susceptible to homogeneous pair¬ 
ing. Tc is then obtained via (the Thouless criterion [30]) 

_ _]_ { i-2rtF(gfc g) _ 2 \ _ r, 

47ras 4V l^k,a „ |fc|2 ) ~ 

The equation of state of the system is determined from 
Eq. (|^ as 




V ^k,a (Cfca) 


where arg(z) is the argument of the c-number z, nsix) = 
l/(e“/^ —1) is the Bose function, and Miuj^, q) is the an¬ 
alytic continuation oiM{iqt z, q) evaluated just above 
the real axis in the z-plane. /i at a given T is then de¬ 
termined from the solution of the equation p{T,fi) = po. 
M{z,q) may have an isolated zero (below the scattering 
threshold a;o(q)) at 2 ; = oJhiq) along the real axis; this sig¬ 
nals the presence of a bosonic bound state of a fermion 
pair with center of mass momentum q. It is useful to 
rewrite the equation of state by explicitly identifying the 
contributions to p 


P{T, p) = Pf{T, p) + Pb{T, p) + pc{T, p) 


(5) 


where pf{T,p) = yJ2k,a'^F i^^a) is the fermion con¬ 
tribution, Pb{T,p) = is the 

contribution from the bosonic poles, and pc(T,p) = 

-V nB[uj) -^ is the contri¬ 

bution from the scattering continuum manifested as a 
branch cut of M{z,q) along the real z-axis. 
Inadequacy of the Gaussian Theory: The Gaussian 
theory, notably successful [28] in the description of the 
interacting Fermi gas in free vacuum, has rather pecu¬ 
liar features in the presence of RSOC (non-Abelian gauge 
field) of the type A = (Ar, A,., Ap) where A^. > Ap. While 



Figure 1. (Color online) Dependence of chemical poten¬ 
tial (/i) on RSOC strength: The Gaussian theory has two 
distinct solutions for p, called the free vacuum branch (FB) 
(higher p, solid line), and the rashbon branch (RB) (lower p, 
dashed line). RB always has the lower free energy within the 
Gaussian theory. The solid line with dots indicates p in the 
approximate non-Gaussian (ANG) theory developed in this 
paper. Non-Gaussian effects eliminate RB for X < kp- 


we focus on the spherical gauge field, our discussion will 
be applicable to all such gauge fields. 

Fig. shows the dependence of /r on A at a fixed low T 
(T < Tp) and negative as- The remarkable feature is that 
the Gaussian theory has two solutions for p for a given set 
of parameters. For X kp, one of the solutions, called 
the “free vacuum branch” (FB) (see Fig. [^, is smoothly 
connected to that of the free vacuum found in previous 
works. pSl The other solution always has p < 0, even 
for X kp. For large A, p in this branch is determined 
by the rashbon dispersion, and hence called the rashbon 
branch (RB). Guriously p along RB, which always has the 
lower free energy, approaches 0“ as A —>■ 0. This suggests 
that the equilibrium state of the Gaussian theory with 
RSOC is not continuously connected to the free vacuum 
in the limit of A —>■ 0! (see Fig. [^) 

The physics of RB at small A can be traced to the 
contribution pb from the bound bosonic states to the 
total density (Eq. (§)• Fig. I shows the dispersion of 
such bosons as a function of their momentum q. Key 
points to be noted, whenever p < 0, are (i) even for neg¬ 
ative scattering lengths, there are bound bosonic states 
in RB whenever |q| < qo(= ■^)) while they cease to 
exist at larger |q|. (ii) the binding energy of the bosons 
{Eb{q) = a;o(q)—Wb(q)), even though significant for small 
|q|, is vanishingly small in the range ® < |q| < qg- This 
physics is quite similar to what is found in the two body 
problem. |21j Such a bosonic dispersion can therefore ac¬ 
commodate a large number of particles forcing p to be 
self-consistently negative. Although E;,(q = 0) —0 for 
vanishingly small A this phenomenon persists, resulting 
in RB not being smoothly connected to the free vacuum. 
Note also that FB does not have any contribution from 
Pb, since in this regime there is no bosonic bound state 
for p > 0. 
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Figure 2. (Color online) Energy dispersion of the bonnd 
boson (fermion-pair) in the Gaussian theory: The 

dashed black line is the scattering threshold uJo{q) as a func¬ 
tion of \q\ (momentum of the pair). 



Figure 3. (Color online) Regime where non-Ganssian ef¬ 
fects are ineluctable: At a given temperature T and scat¬ 
tering length fls, non-Gaussian effects play a crucial role when 
A < Xng{T,as). Inset shows the parameter b (see text) that 
characterizes the non-Gaussian effects at a fixed temperature. 


As is evident, the fully formed bound states in the 
range ^ ^ |<7| < 9o with a vanishing binding energy can 
easily be destabilized. In particular, the quartic term in 
Eq.§ with the coupling K describes the interactions be¬ 
tween the pairing fluctuations r]{q). A natural question 
is whether these weakly bound states are stable when the 
interactions between the bosonic fluctuations are taken 
into account. We address this question by quantifying 
the strength of these beyond-Gaussian effects by the pa¬ 
rameter b (proportional to K{qi, 92 ; 93 , 94 ) in the limit of 
zero momentum |251 [27] ) obtained as 


_ J_ f 1 —2 rtf (Cfc a) 

4V 2^k,OL 


2 nF{^ka)0- — '<^F{^ka)) 


( 6 ) 


When /r is large and negative, as in a “boson dominated” 
state where the most prominent contribution arises from 
Pi,, b Ki 327 rV 2 (-|)) 5/2 • physical meaning of b can be 
made evident by noting that b ~ when A = 0 and 

the scattering length is small positive (free vacuum BEG 
side). Here obb is the scattering length of two bosons 
(bound fermion pairs) and is proportional to a..|22j Fur¬ 
thermore, for any Og, as A —)■ 00 , b —>■ A“^, which can 
be immediately identified with where u/jb is the 

rashbon-rashbon scattering length. p5] Therefore, b^^^, is 
a length scale that characterizes the interactions among 
the pairing fields. Interestingly, this parameter is nonzero 
in the limit of A —>■ O'*" and grows with increasing A at¬ 
taining a peak when A « kp{see inset of Fig. |^, subse¬ 
quently possessing the just discussed asymptotic behav¬ 
ior at large A. 

The effects of b on the weakly bound states can now 
be estimated in a physical manner. The lowest order ef¬ 
fect of b would be to shift the energy of the bound state 
via a Hartree shift, i.e., uJb{q) —>■ ujb{q) + p) 

where k is a dimensionless number of order unity [31]. 
Clearly, the bound bosonic state will be unstable if the 
shift takes it into the scattering continuum, i.e., a neces¬ 
sary condition for the stability of the bound state is that 


ujb{q) + nb^/^pb{T, p) < uJo{q). Thus RB can be stable 
only if 

Eb{q = 0)>Rb^/^PbiT,p). (7) 

Using this criterion we obtain the regime (see Fig. 
where RB is eliminated by non-Gaussian effects i.e., 
A < Xng{T, Os). For a given Og, we find that Xng increases 
with decreasing temperature |32j . These estimates pro¬ 
vide a lower bound of Xng, which results in Xng ~ kp for 
temperatures T <Tp. Thus beyond-Gaussian effects are 
crucial in the most interesting regime of parameters. 
Approximate Non-Gaussian Theory: Having hrmly 
established that even a qualitatively correct description 
of spin-orbit coupled Fermi gases necessarily requires a 
beyond-Gaussian theory, we propose and discuss one such 
theory. A key desideratum of such a theory is the elim¬ 
ination of RB for A ^ Xng, and a smooth evolution (at 
given T, Og) from the free vacuum state at vanishing A 
to the rashbon gas at large A. The implementation of 
such a theory is a formidable challenge, even as we note 
that Gaussian theory itself requires considerable calcula- 
tional effort [33|. Faced with this reality, we develop an 
approximate non-Gaussian (ANG) theory by a suitable 
modification to the equation of state (Eq. Q) that only 
entails the same calculational complexity as the Gaus¬ 
sian theory. The approximation follows the physical ar¬ 
gument that the non-Gaussian term b shifts ujbiq) to 
^b{q) + nb^/^Pb^ where p^^ is the bound boson con¬ 
tribution calculated within the Gaussian approximation. 
Only those bosonic states that remain below the scat¬ 
tering continuum after this energy shift, i.e., the bosonic 
states for all |< 7 | < qb obtained by Eb{qb) = nb^^^pf^, are 
stable to non-Gaussian effects. These arguments pro¬ 
vide for the approximation to the contribution of the 
bosonic bound pairs to the equation of state (Eq. @) 
as p^^^{T,p) = -yT,\g\<g, riBiuJbiq))^^^- 

Fig. [T] shows the results of this approximate non- 
Gaussian theory (see the curve marked ANG). The ap- 
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Figure 4. (Color online) Phase diagram: Dependence of 
superfluid Tc on aa and A (ANG theory). Crossover occurs 
in regime enclosed by the dashed lines. The dotted line indi¬ 
cates a candidate path that can be traced out in a cold atoms 
experiment by varying the density at a fixed scattering length 
and RSOC. 


proximate theory does indeed possess the key features 
desired, notably the elimination of RB for A <C /cf- In 
this regime, the ANG chemical potential smoothly con¬ 
nects to free vacuum value. Furthermore, when A 3> /cf, 
the ANG recovers the rashbon gas. In the ANG theory, 
RB appears only after a particular value of A which de¬ 
pends on T and Og at which the solution switches from 
FB to RB. This evolution should be smooth in a detailed 
theory which also includes non-Gaussian effects in the 
free vacuum. 

Phase Diagram: We now use the ANG theory to obtain 
the phase diagram shown in Fig. For a small positive 
scattering length, which obtains a BEG of fermion pairs 
in free vacuum, increasing RSOG engenders a smooth 
crossover to the rashbon-BEG with the Tc gradually 
changing from that set by the free vacuum boson mass 
(twice the fermion mass) to that set by the rashbon mass. 
At the resonant scattering length, the Tc again evolves 
from that of the free vacuum unitary Fermi gas, to that of 
the rashbon-BEG. The scenario for a small negative scat¬ 
tering length is significantly different as discussed below. 

One of the key aspects of the phase diagram Fig. 
shown in detail in Fig. is the large enhancement of 
Tc for a system with a weak attractive interaction. For 
example, for = 2, the Tc is enhanced from 0.02Tf 
at A = 0 to about O.ITf when A « kp. Further, there 
is a regime of A where Tc decreases. Beyond this, Tc is 
determined by two-body physics as shown by dashed-dot 
lines in Fig. H Fig. [5] also shows the mean field Tc which 
makes fluctuation effects evident. For example, the 
from the ANG theory is about 85% of the mean field Tc 



Figure 5. (Color online) Non-monotonic dependence of 
Tc on A for weak attraction: Points indicate the Tc cal¬ 
culated from our ANG theory (lines through the points are 
a guide to the eye). Dashed lines - Tc from mean field the¬ 
ory. Dashed-dot lines - Tc estimated from the condensation 
of the bound-fermion pairs. Thin horizontal line - Tc of the 
rashbon-BEG. 

for —-r^ = 2.0, and it is reduced to about 60% of the 
mean field value when --r-^ — = 0.75. The enhancement 

Kfras 

of Tc for weak attractive interactions indicated by our 
ANG theory is a remarkable feature of spin orbit coupled 
systems (non-Abelian gSMge fields). Finally, for any given 
Os (including Og < 0), the Tc at large A (A > (kp, l/|ag|)) 
is independent of Og, determined only by the rashbon 
mass. 

As shown in Fig.|^ there is a much bigger regime of pa¬ 
rameters (with weak interactions and RSOG) over which 
the crossover from a BCS like ground state to a rashbon- 
BEG occurs. The central point is that the superfluid 
with high Tc occurs in this crossover regime. Indeed, it 
will be interesting to mimic this crucial finding in mate¬ 
rial systems to provide routes to making superconductors 
with high transition temperatures. On a different token, 
this physics can be uncovered in a cold atoms experi¬ 
ment at fixed negative Og and RSOG, by working with 
different trap centre densities, tracing out a path akin 
to the dotted line shown in Fig. Another interesting 
point to note is that the enhanced binding induced by 
the RSOG will result in significant pseudogap features |M] 
which could be observed even at higher temperatures. 

In summary, we have shown the crucial role of beyond- 
Gaussian effects in spin orbit coupled Fermi gases. We 
have developed a simple theory that incorporates the 
beyond-Gaussian effects in an approximate fashion. Us¬ 
ing this theory we obtain the phase diagram of the sys¬ 
tem. A key result of our calculation is the demonstra¬ 
tion of the enhancement of the exponentially small su¬ 
perfluid transition temperature with weak attraction to 
values comparable to Fermi temperature. This impor¬ 
tant point provides clues to producing superconductors 
with high transition temperatures. Our approximate 
non-Gaussian theory uncovers the rich physics in spin- 
orbit coupled gases providing motivation for further de- 
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tailed theoretical considerations. Promising routes to 
treat beyond-Gaussian effects include the Gq —G or G—G 
schemes. [551 [5S] 
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relative momentum of a two-particle state with particles 
having helicities a and k) is the weight of such 

a state in the singlet sector. The third term in Eq. ( |A1[ ) 
contains external pairing source fields F(q), 


= V ? E n<l)S*{q) + F^{q)S{q). (A4) 

Appendix A: Fluctuation Theory Formulation ' „ 


Many body physics is studied by introduc¬ 
ing four component Nambu fields = 

{c\(k) c+(—fc) c*_{k) c_(—fc)), where k denotes the 
four-vector {ikn,k), ikn being a Fermi-Matsubara fre¬ 
quency. The action of a system of interacting fermions 
is made up of three pieces: 


This term anticipates a pairing instability in the system, 
and is added solely to aid the calculation of the pairing 
susceptibility (most of the formulae, therefore, will have 
F = 0). 

We now perform a Hubbard-Stratanovich transforma¬ 
tion on Sv by introducing pairing fields r]{q), 


-l--l-(-^1) 

The term Sq describes the kinetic energy including the 
RSOC induced by the non-Abelian gauge field, 

M'l'] = ^ E •^*(fc)(-Go'(fc, k'))F{F). (A2) 


where GQ^{k,F) = Diag(f/c„ - + -I- ^k + ,ikn - 

+ Cfe-) 4 .fc' with ^ka = £ka- and fj, is the 
chemical potential. The second term, S^,, in the ac¬ 
tion (Eq. (AD) describes the contact attraction among 
fermions as 


= (A3) 

9 

where T and V are respectively the temperature and vol¬ 
ume of the system, v is the bare interaction parame¬ 
ter. This last quantity is traded for the s—wave scatter¬ 
ing length Us through regularization as = ^ -|- A, 
where A = ^ J^k W denoting the ultraviolet cutoff. The 
quantity 5'*(g) = k) {^ + k) - k) 

stands for the singlet pair density in Matsubara-Fourier 
space, q = {iqe,q), where iq^ is the Bose-Matsubara fre¬ 
quency, q is the center of mass momentum and k is the 


S[F, v,F]=Y, F*{k){-G-\k, ^*(9) V{q), 

k,k' Q 

(A5) 

where 

G-i(fc, k') = Go i(fc, k') - 7(fc, k'), (A 6 ) 


'y{k,k') 


0 7++(/c,fc') 0 7+_(fc,fc') 

7++(fc,fc') 0 7+_(fc, fc') 0 

0 7 ^(fc, fc') 0 7 {k,k') 

7 y{k,k') 0 7 {k,k') 0 

(A7) 


with 7 c,/ 3 (fc,fc') = T,g'yiq)Aci3{q,k- ^)6q^k-k', 

%l3{k,k') = T,qy*{-q)A*fi^{-q,k- ^)Sg^k-k', and 

liq) = ’qiq) + F{q). The action is now quadratic in 
fermionic fields which can be integrated to yield 


=-lndet[-G - V r?*(g) r?(g). (AS) 

V ^^ 

Q 


This is Eq. 0 in the main text. 



